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SEMIANNUAL  REPORT  SUMMARY 


The  four  projects  being  carried  out  in  the  area  of  Laser  Processes 
under  this  contract  are  summarized  below.  More  detailed  discussions  are 
given  in  Sections  II  through  V of  this  report. 

(1)  Metal  Vapor-Rare  Gas  Discharges. 

Several  metal  vapor-rare  gas  molecules  have  been  proposed  as 
the  active  media  for  high  efficiency,  high  power  exclmer  lasers 
operating  at  near  visible  wavelengths.  Our  previous  research  in  this  area 
had  been  directed  primarily  toward  the  measurement  of  the  stimulated 
emission  and  absorption  coefficients  for  some  of  these  molecules.  The  present 
work  is  a more  detailed  evaluation  of  the  potential  for  laser  use  of 
one  of  the  more  promising  of  these  molecules,  i.e.,  the  NaXe  molecule. 

In  order  to  take  advantage  of  the  high  quantum  efficiency  inherent  in 
these  excimer  systems,  i.e.,  greater  than  80%,  it  is  desirable  to  excite 
the  associated  atomic  levels  as  efficiently  as  possible.  We  believe 
that  electric  discharge  excitation  with  its  low  mean  electron  energies 
offers  the  best  possibility  for  efficient  excitation  and  so  have  initiated 
the  present  investigation  of  electric  discharges  in  mixtures  of  Na  and  Xe 
at  high  Xe  densities.  The  effort  covered  by  this  report  includes 
a)  the  development  and  application  of  models  of  the  discharge  and  its 
radiative  properties  and  b)  the  construction  and  initial  operation  of  a 
small  volume  discharge  equipped  with  apparatus  for  the  measurement  of 
the  electrical  and  radiative  properties  of  the  discharge.  The  principal 
specific  conclusions  drawn  from  the  modeling  calculations  during  this 
contract  period  are  that  presence  of  significant  excited  state  populations 
and/or  of  a high  radiative  flux  level  drastically  changes  the  electron  energy 
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balance and  the  degree  of  ionization,  that  the  earlier  indication  of  a large 
reduction  Na2  density  by  electron  collisional  dissociation  is  correct, 
and  that  the  model  is  incomplete  and  will  have  to  be  extended  to  include 
at  least  several  more  excited  states  of  the  sodium  atom.  The  general 
conclusion  is  that  discharge  excited  metal  vapor-rare  gas  systems 
continue  to  appear  to  be  viable  candidates  as  high  efficiency, 
high  power  lasers  for  near  visible  wavelengths. 

(2)  Stability  of  Discharges  in  Weakly  Ionized  Gases. 

The  technical  problem  being  investigated  is  the  optimization  of 
the  use  of  electrical  energy  for  the  production  of  excited  molecules 
in  electrically  excited  gas  lasers.  In  particular,  the  efficient  use 
of  electrical  excitation  in  gas  lasers  requires  that  the  discharge  remain 
diffuse  as  the  density  of  the  excited  atoms  and  molecules  is  raised  as  high  as 
possible,  i.e.,  that  the  discharge  not  form  a constricted  channel  or  arc. 

There  is  evidence  that  arcing  in  lasers  originates  in  or  near  the  region 

of  very  high  energy  input  at  the  cathode,  known  as  the  cathode  fall.  During  this 

report  period  we  have  essentially  completed  our  calculations  of  the  time 

development  of  the  cathode  fall  region  of  a high  oressure  discharge  and 

have  initiated  an  analysis  of  published  experimental  data  based  on  the 

scaling  laws  derived  from  the  theory. 

(3)  Electron  Excitation  of  Molecular  Metastables. 

The  electron  excitation  of  molecular  metastables  in  an  electrical  discharge 
can  be  an  efficient  means  of  producing  high  excited  state  densities  in 
gas  lasers.  During  this  report  period  we  have  extended  the  measurements  of 
electron  excitation  rate  coefficients  for  the  (^(b^X*)  metastable  state  to 
excitation  in  mixtures  of  O2  and  N2  and  have  initiated  measurements  of  electron 
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rate  coefficients  for  the  aV  state  of  N2  in  pure  N2.  The  O2  (b1!) 
measurements  in  mixtures  confirm  our  results  in  pure  i.e., 

the  measured  total  excitation  rate  coefficients  at  mean  electron  energies 
of  about  2 eV  are  as  much  as  a factor  of  three  larger  than  expected 
on  the  basis  of  previously  available  data. 

(&)  Scattering  and  Transport  of  Resonance  Radiation  in  Gases. 

The  technical  problem  considered  here  is  the  measurement  of  the  transport 
and  scattering  of  resonance  radiation  emitted  by  metal  vapors  so  as  to  obtain 
the  rates  of  radiative  and  non-radiative  energy  loss  by  excited 
atoms  in  excimer  lasers,  etc.  The  measurements  of  the  spectral  intensity 
of  scattered  radiation  under  white  light  illumination  have  been  completed 
and  initial  comparisons  have  been  made  with  theory.  These  measurements 
confirm  the  loss  of  excited  Na  atoms  by  excitation  transfer  to  ^2  molecules. 
The  integrated  scattered  intensities  confirm  the  validity  of  theoretical  model 
at  low  (<10  atom/cm  ) and  at  high  (>10i  atom/ cm  ).  The  discrepancies  at 
intermediate  Na  densities  may  be  caused  by  the  approximate  manner  in  which 
the  theory  treats  coherently  scattered  radiation.  Also,  during  this  report 
period  experiments  were  initiated  in  which  a dye  laser  will  be  used  to  excite 
the  Na  resonance  state  near  line  center  so  as  to  emphasize  the  effects  of  the 
non-radiative  excitation  transfer  process. 

METAL  VAPOR-RARE  GAS  DISCHARGES 

Drs.  H.  Rothwell,  R.  Shuker,  and  A.  Gallagher  and  Mr.  D.  Leep. 

The  objective  of  this  project  is  to  determine  whether  metal  vapor-rare 
gas  discharges  can  be  operated  so  as  to  produce  the  degree  of  metal  vapor-rare 
gas  excimer  excitation  required  for  high  power  visible  laser  opeationn. 

The  project  includes  experimental  measurements  and  theoretical  models  of  the 
properties  of  Na-Xe  discharges  at  rare  gas  pressures  up  to  about  10  atm. 


. v 


# A* 


-6- 


The  Na-Xe  mixture  is  chosen  because  we  have  more  of  the  data  needed 

for  a meaningful  model^  and  because  this  mixture  might  yield  a high  power 

visible  excimer  laser.  In  particular,  our  previous  calculations  indicate 

2 

that  potentially  ruinous  effects,  such  as  absorption  by  ground  state  Na^ , 

are  not  important."^  The  following  discussion  summarizes  our  modeling 

efforts  up  to  the  date  of  this  report.  The  portion  of  the  work 

carried  out  under  this  contract  includes  a)  the  development  and 

application  of  procedures  for  the  iteration  of  the  solutions  for  the 

electron  energy  distribution  and  for  the  populations  of  the  excited  states, 

ions,  and  electrons;  b)  the  inclusion  of  the  effects  of  stimulated  emission 

and  of  the  intercavity  flux  in  the  model;  and  c)  the  effects  of  the  addition 

of  a third  excited  state  on  the  discharge  gain,  etc. 

The  A-X  bands  of  alkali-noble  gas  excimers  and  alkali-dimers  have  been 

proposed  for  high-power,  efficient  lasers  and  gain  cells  in  the  650-1300  nm 
2,3 

wavelength  region.  The  excited  states  of  these  excimers  are  weakly 

4 

bound  compared  to  those  of  the  rare  gas  molecules,  but  net  gain  was 
predicted  with  a fev  percent  excited  alkali  atoms  due  to  a strongly 
repulsive  ground  state.  The  alkali  dimers  also  yielded  net  gain  on  the 
A-X  bands  at  a few  percent  excitation,  due  to  displaced  A and  X state 
potentials  which  yield  emission  centered  at  longer  wavelengths  than 
absorption. 

A typical  gas  mixture  considered  here  is  10^  - 10^  atoms/cm^  of 

20  3 

alkali,  with  a few  percent  bound  as  dimers,  and  about  10  atoms/cm  of 
noble  gas  to  act  as  a heat  sink  and  a constituent  of  the  excimer  molecule. 
Since  the  excimer  and  dimer  bands  are  associated  with  the  first  excited 
alkali  state,  while  the  noble  gas  excited  states  are  at  much  higher 


energies,  electric  discharge  excitation  is  an  attractive  possibility. 
Simple  models  of  a laser  utilizing  K-Xe  have  been  given  by  Palmer^ 
using  Maxwellian  electron-energy  distributions  f(u)  and  a limited 
sot  of  population  rate  equations.  Schlie^  has  shown  that  at  low 
fractional  ionization  the  rapidly  decreasing  non-Maxwellian  tail  of  f(u) 
allows  efficient  excitation  of  the  lasing  alkali  state  rather  than  the 
noble  gas.  The  present  work  has  modeled  this  laser  by  solving  self- 
consistently  the  Boltzmann  equation  for  f (u)  and  an  extensive  set  of  rate 


equations  for  the  various  excited  and  ionized  species  densities.  We  assume 
a homogeneous  self-sustained  discharge  at  fixed  E/N.  The  questions  of 


instabilities  or  constrictions  are  being  studied  experimentally  in  this 
7 

laboratory. 


We  have  modeled  a Na-Xe  mixture  since  the  Na£  and  NaXe  A-X  bands  are 
2 8 9 

fairly  well  understood,  ’ * many  electron-Na  cross  sections  are 
available, 10, and  the  population  inversion  obtained  with  Na-Xe  is 


higher  than  with  Na  and  other  noble  gases.  We  report  here  calculations 

for  the  case  of  [Na]  = 3.  5*  10^  cm  [Xe]  = 2.7*10^  cm  T = 760  K, 

—18  2 

and  E/N  = 8 * 10  V cm  , where  [ ] indicates  species  density,  T is  the  gas 


temperature,  E is  the  applied  electric  field,  and  N is  the  total  gas 


density.  Figure  1 shows  the  Na,  Na^,  and  NaXe  levels  included  in  the 

model,  as  well  as  the  most  important  processes.  Due  to  close  coupling 

by  electron  collisions  the  Na  4S  and  3D  levels  are  combined  into  a "4S" 

level,  and  the  Na^(A  E ) and  Na2(a  II)  are  combined  into  an  A'  level 
1 3 2 

composed  of  1/7 (A  X)  state  and  6/7(a  IT)  state.  Electron  momentum  transfer, 

excitation,  deexcitation,  ionization,  and  dissociative  recombination  are 

11  —3 

included.  The  calculation  starts  with  an  electron  density  of  ng  = 10  cm 
and  typically  reaches  steady  state  in  a few  microseconds.  At  steady-state 
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the  rate  of  ionization  of  the  excited  states  Na(3P)  and  Na("4S")  is  60  and 

140  times  larger,  respectively,  than  that  from  ground  state  sodium. 

—16  2 

Excitation  and  ionization  of  Xe  plays  a minor  role  at  E/N  <10  V cm  . 

Rapid  association  of  Na+  to  form  molecular  ions  yields  NaXe+  as  the  dominant 
ion.  Dissociative  recombination  balances,  at  steady  state,  the  production 

of  ions  and  electrons.  The  rate  coefficients  for  dissociative  recombin- 

+ + - 1/2  -7  3 -1 

ation  of  NaXe  and  Na^  , are  assumed  to  be  (0.04/u)  *10  cm  sec  , 

where  u is  the  average  electron  energy  in  eV.  Increasing  this  value  raises 
the  E/N  required  for  a given  electron  density  and  gain.  Dissociative  recombin- 
ation is  assumed  to  populate  equally  the  Na  3P  and  "4S"  states. 

The  ratio  [NaXe (A) ] / [Na(3P) ] remains  very  near  the  equilibrium  fraction 
(0.2  for  this  [Xe]  and  gas  temperature)  since  the  weakly  bound  A state  is 
rapidly  formed  and  dissociated  by  Xe  collisions.  Thus,  the  gain  in  the 
excimer  band  is  fixed  by  [Na(3P)J,  which  is  in  turn  fixed  predominately  by 
the  ratio  of  3S  •>  3P  excitation  to  deexcitation  by  electrons.  Radiative 
decay  rates  include  radiation  trapping  of  the  Na  resonance  lines. 

In  the  absence  of  excited  Na,  and  thereby  deexcitation  collisions,  the 
calculated  f(u)  agree  with  those  of  Schlie^  for  the  same  conditions.  However, 
including  collisional  deexcitation  of  the  excited  Na  raises  f(u)  at  energies 
above  the  first  excitation  energy  of  sodium  at  2.1  eV.  In  this  case 
deexcitation  shifts  more  electrons  to  u > 2.1  eV  than  are  removed  by 
excitation  and  ionization  of  excited  states.  This  increases  the  3S  -*■  3P 
excitation  and  excited-state  ionization  rate  coefficients. 

It  was  noted  in  Ref.  2 that  absorption  by  the  Na2>  X to  A band  could 
cancel  the  excimer  gain  if  [Na2(X)]  were  not  depleted  by  thermal  or  other 
mechanisms.  One  of  the  principal  conclusions  of  the  present  study  is  that 
[Na2(X)]  is  greatly  depleted  by  electron  collisional  dissociation  at  the 
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steady-state  ng  appropriate  for  high-power  operation.  This  dissociation 

process  occurs  via  electron  collision  excitation  from  X3 1 to  the  repulsive 

g 

3 

x E state,  which  dissociates  to  two  Na(3S)  atoms.  The  identical  process 
12  13  - 

has  been  measured  and  calculated  for  dissociation.  The  estimated 


Na_ 


2 cross  section  is  obtained  by  multiplying  the  measured  cross  section 
2 

by  p and  dividing  its  energy  scale  by  p,  where  p is  the  ratio  of  the 

excitation  thresholds  for  and  Na2  • This  form  of  scaling  has  been  found 

accurate  in  many  cases. 34  3^  A similar  process  with  twice  the  cross  section 

is  assumed  to  deplete  the  ^2^')  state.  Depletion  resulting  from  gas 

heating  and  electron-induced  vibrational  excitation  are  not  included. 

Thus  the  actual  depletion  should  exceed  that  calculated 

For  the  above  total  Na  and  Xe  densities  and  E/N,  the  rate  equations 

without  inclusion  of  stimulated  emission  yield  steady-state  values  of 

n = 3 x 1015  cm"3,  [Na(3P>]  = 3 x 1015  cm"3,  [Na("4S")]  = 4 x 1034  cm"3, 

6 

[Na  (X)]  = 9 x 1013  cm"3,  faa2(A)]  = 1/7  x [Na2 (A' ) ] = 4 x 1012  cm"3, 

[Na+]  = IQ34  cm  3,  [Na2+]  = 1034  cm  3,  and  [NaXe  ] = 3 x 10  3 cm  . The 
resulting  steady-state  stimulated  emission  g and  absorption  coefficients 


are  given  in  Fig.  2.  These  curves  are  based  on  the  normalized  and  g^ 
coefficients  given  in  Refs.  2 and  8,  except  that  the  radiative  transition 
probability  for  the  Na.  A-X  band  is  twice32  that  used  in  Ref.  2.  The  Na„(X) 

V ^ ^ 

absorption  shown  in  Fig.  2 is  smaller  by  a factor  of  20  than  that  predicted 

-1  -3  -1 

by  thermal  equilibrium.  A net  gain  of  10  - 10  cm  is  indicated  for  tfte 

-3  -1 

excimer  band.  The  Na2  A-X  band  has  only  ~10  cm  gain  in  this  figure, 

but  can  be  considerably  larger  at  lower  E/N  and  n^,  where  less  depletion 

occurs.  The  excimer  transition  is  homogeneously  broadened  on  a 10  33>  sec 

18 

time  scale,  as  determined  by  measured  excimer  dissociation  rates.  Thus, 
the  laser  could  operate  at  any  wavelength  from  680-900  nm. 


STIMULATED  EMISSION  AND 
ABSORPTION  COEFFICIENTS  (cm-1) 
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X ( nm ) 


Fig.  2.  Calculated  stimulated  emission  and  absorption  coefficients  for 
[Na]  = 3.5X1016  cm”3,  [Xe]  = 2.7*1020  cm'3,  T = 760  K,  and 
E/N  = 8x10”  V-cm  . As  discussed  in  the  text  the  result  of  this 
preliminary  model  must  be  used  with  caution  pending  a thorough 
investigation  of  the  effects  of  the  addition  of  higher  excited  states. 
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The  calculations  show  that  the  homogeneous  discharge  has  a positive 
volt-ampere  characteristic  with  a resistivity  of  about  20  fl-cm.  This 
behavior  is  the  result  of  a balance  between  ionization  via  excited  Na 
atoms,  which  increases  approximately  linearly  with  electron  density,  and  the  loss  of 
electrons  and  ions,  which  is  proportional  to  the  square  of  tne  electron 
density. 

The  rate  equations  have  also  been  solved  with  the  inclusion  of 

stimulated  emission  to  evaluate  the  potential  laser  power  and  efficiency. 

As  the  radiative  flux  is  increased  E/N  must  be  raised  so  as  to  maintain 

sufficient  excited  state  and  electron  density.  Thus,  Fig.  3 shows  the 

dependence  of  the  net  gain  at  X = 700  nm  on  the  cavity  power  at  steady- 

-17  2 

state  for  E/N  = 1 x 10  V cm  . It  also  shows  the  corresponding  radiative 
power  output  and  electrical  power  input.  The  solid  curves  of  Fig.  3 show 
calculations  for  the  full  value  of  the  estimated  cross  sections  for 
electron  impact  dissociation  of  the  Na2(X)  and  Na2(A')  states,  while  the 
dashed  curves  are  for  one  tenth  of  these  cross  sections.  The  presence  of 
a high  cavity  flux  causes  bleaching  of  the  Na2(A-X)  band  reducing  the 
absorption. 

As  an  example  of  laser  characteristics  predicted  from  the  calculations 

of  Fig.  3,  consider  the  conditions  for  the  case  of  the  full  Na2  dissociation 

cross  section  and  the  maximum  radiant  power  output.  Using  the  spatially 

19 

uniform  flux  approximation  for  a laser  oscillator  with  a single 

transmitting  mirror  we  find  that  if  the  laser  were  50  cm  long  and  a 

50%  transmission  output  mirror  were  used,  the  output  power  would  be  about 
2 

1.5  MW/cm  and  the  overall  efficiency  would  be  about 


GAIN  COEFFICIENT  (cm- 


CAVITY  INTENSITY  (MW/cm2) 

Fig.  3.  Calculated  gain  and  power  output  vs.  cavity  flux  at  700  nm 

for  Na-Xe  discharge.  See  text  for  rp'  arks  on  the  validity  of 
this  preliminary  model. 


POWER  (MW/Jt) 
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10%.  The  intercavity  flux  would  be  about  5 MW/cm  and  the  input  power 
would  be  30  MW/S,.  The  specific  heat  of  the  Xe  results  in  a gas  temperature 
rise  of  100  K/usec,  so  that  the  laser  pulse  might  terminate  after  a few 

microseconds  of  operation.  These  conditions  have  not  yet  been  optimized 

* 

with  respect  to  the  output  coupling,  the  [Na]/[Xe]  ratio,  the  [Xe]  density 
or  E/N. 


This  initial  modeling  of  electric  discharge  excitation  of  the  Na-Xe 

excimer  indicates  that  this  system  could  yield  an  efficient,  tunable, 

high-power  laser  or  gain  cell  operating  near  700  nm.  Obviously  a number 

of  practical  problems  remain  to  be  addressed,  e.g.,  the  optical  quality 

of  the  medium  and  the  practical  operation  of  Na-Xe  discharges  at  the  high  E/N 

of  these  models.  The  principal  uncertainty  in  the  model  appears  to  be 

whether  the  use  of  two  or  three  excited  states  adequately  represents  the  effects 

20 

of  multistep  ionization.  Results  obtained  subsequent  to  this  report  period 
show  that  a correct  model  will  somehow  include  the  effects  of  a large  number 
of  excited  states. 

The  experimental  portion  of  this  effort  included  the  design  and  construction 

of  an  improved  high  temperature  discharge  cell,  the  calibration  of  the  optical 

system  for  absolute  intensity  measurements,  and  two  preliminary  spectral  scans 

at  relatively  low  discharge  current  densities.  Figure  4 shows  the  spectral 

distribution  of  radiation  emitted  at  [Na]  = 7.5x10^  cm~^,  [Xe]  = 4.6x10'^  cm  ^ 

2 

and  ^ 100  A/cm  . The  dashed  curves  show  the  contributions  of  the  NaXe  (A+X) 
band  and  the  Na£  (A+X)  band.  Note  that  at  these  high  current  densities 
the  Na2  band  intensity  is  much  lower  than  expected  for  an  equilibrium  ratio  of 
Na2  to  Na,  i.e.,  Na2  depletion  is  large  as  predicted  by  our  model  for 

9 

low  current  densities.  During  the  next  report  period  the  absolute  intensities 
of  the  various  bands  and  lines  will  be  used  to  calculate  Na  and  NaXe  excited 
state  densities  for  comparison  with  discharge  models  and  to  predict  discharge 
characteristics  at  higher  power  densities. 


Relative  spontaneous  emission  per  unit  wave  number. 
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Wavelength,  nm 


Fig.  4.  Experimental  spectral  intensity  distribution  from  Na-Xe  discharge 
for  [Na]  = 7.5*1015  cm"3,  [Xe]  = 4.6*1019  cm-3  and  I % 13.4  A. 

pK 
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III.  STABILITY  OF  DISCHARGES  IN  WEAKLY  IONIZED  GASES 
Drs.  H.-C.  Chen  and  A.  V.  Phelps 

The  objective  of  this  theoretical  investigation  of  the  growth  of 
instabilities  in  weakly  ionized  gas  discharge  is  to  develop  quantitative 
techniques  for  the  prediction  of  conditions  for  formation  of  the  arcs 
which  limit  the  power  output  from  high  power  laser  systems.  The  avail- 
ability of  such  techniques  would  greatly  aid  the  optimization  of  high  power 
lasers  currently  under  development  and  would  greatly  improve  the  reliability 
of  the  design  of  scaled-up  versions  of  these  lasers.  During  this  contract 
period  we  have  essentially  completed  our  development  of  techniques  for 
the  detailed  calculation  of  the  growth  of  the  cathode  fall  region  of  a 
glow  discharge.  The  model  described  in  detail  below  uses  one  dimensional 
solutions  of  the  continuity  equations  for  electrons  and  positive  ions  and 
Poisson's  equation.  The  electron  and  ion  transport  coefficients  and  the 
ionization  and  recombination  coefficients  were  based  on  experimental  measure- 
ments for  helium,  which  is  used  as  the  principal  constituent  in  most 
electric  discharge  lasers.  These  coefficients  are  readily  modified  to  allow 
the  modeling  of  gas  mixtures.  Thus  far  we  have  not  modeled  the  cathode  fall 

for  discharges  using  attaching  gases.  These  calculations  represent  improvements 

21  23 

over  previous  analyses  in  that  the  earlier  work  omitted  electron  diffusion  ’ 

22  23 

and/or  presented  only  very  limited  results.  ’ 

A.  Theory 

(1)  Basic  Equations 

The  electron  and  positive  ion  rate  equations,  Poisson's  equation,  and 
experimental  transport  and  rate  coefficients  are  used  for  tracing  the  growth 
In  time  and  space  of  a one-dimensional  helium  discharge  where  the  electric 
field  is  distorted  by  the  space  charge  of  the  electrons  and  ions.  The 
continuity  equation  for  electrons  is 
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-jf-  + 7 * fcgu)  = vn0  - angn+  , (1) 

where  ng  and  n+  are  the  electron  and  positive  ion  densities,  t is  the  time, 
v is  the  ionization  coefficient,  and  a is  the  recombination  coefficient. 

The  electron  velocity  u can  be  expressed  in  terms  of  drift  and  diffusion 
velocity  as 


V(neD) 

u = - pE » (la) 

e 

where  E is  the  electric  field,  D and  p are  the  diffusion  and  mobility 
coefficients  of  the  electrons,  respectively. 

The  continuity  equation  for  positive  ions  is 
3n+ 

— + 7 • (n+u_))  = vng  - anen+  . (2) 

The  ion  velocity  u+  car.  be  expressed  as 

U+  = , (2a) 


where  p+  is  the  mobility  of  positive  ions.  The  diffusion  of  the  ions  has 
been  neglected  because  of  their  heavy  mass. 

Poisson's  equation  describing  the  space-charge-distorted  fields  is 


3E  = _e 
9z  _ e 


(n. 


ne} 


(3) 


o 

where  e is  the  electronic  charge,  is  the  permittivity  of  free  space 
and  z is  the  position  in  the  axial  direction. 

The  coefficients  D,  v,  a,  u+  and  p must  be  taken  from  the  experimental 
measurement  or  electron  transport  coefficient  calculations.  For  atomic 
helium  they  are  expressed  in  terms  of  the  ratio  of  electric  field  to  gas 
density  E/N  and  shown  in  Eq.  (8)  of  Reference  24. 
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An  allowance  has  been  made  for  the  time  required  for  electrons 
emitted  from  the  cathode  to  attain  the  ionization  rate  appropriate  to 
the  local  E/N  value  by  requiring  the  electrons  to  pass  through 
a potential  of  40  V before  ionization  begins.  See  Druyvesteyn 
and  Penning. ^ 

(2)  Boundary  Conditions 

At  the  cathode,  secondary  electrons  are  produced  by  the  bombardment 
of  ions  and  photons 


r (o)  = y_j\  (0)  + y r . _ . 

e'  '+  + ' p photon’ 


(4) 


where  represents  the  flux  of  various  particles  i.  y+  and  Yp  are  the  secondary 

coefficients  representing  the  probability  of  the  production  of  secondary 

electrons  at  the  cathode  per  incident  positive  ion  or  photon.  The  effect 

of  photons  is  to  enhance  the  production  of  secondary  electrons 

and  speed  up  the  formation  of  the  cathode  fall  at  the  early  times  in  discharge. 

At  the  anode,  we  have  zero  density  for  both  electrons  and  ions 

ne(d)  = n+(d)  = 0 . (5) 

The  voltage  across  the  discharge  gap  is  V so  that 
fd 


'0 


Edz  = V 


(6) 


Equation  (1)  is  a second  order  partial  differential  equation  with  two 
boundary  conditions.  Equations  (2)  and  (3)  are  first  order  partial  differential 
equations  with  one  boundary  condition.  Therefore,  we  have  a complete  set 
of  three  partical  differential  equations  that  can  be  solved  numerically 
for  the  three  unknowns  n , n and  E.  A Crank-Nicholson  finite  difference 

C T* 

26 

scheme  which  is  implicit  both  in  time  and  space  has  been  used  and  a 
predictor-corrector  step  has  been  added  to  take  into  account  the  nonlinear 
coefficients  in  the  equations.  Use  of  nonuniform  spatial  step-sizes 
has  been  found  very  effective  in  studying  the  time  development  of  the 
boundary  layer  in  the  cathode  fall  region. 


f A 
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B.  Results 

As  In  a double  discharge  laser,  an  Initial  electron  and  ion  density 

10  -3 

which  varies  smoothly  with  z and  has  a peak  value  of  3 » 10  cm  is 

assumed.  This  initial  ionization  is  shown  in  Fig.  5.  Calculations  are 

19  -3 

shown  for  a gap  length  of  2 cm,  a helium  density  N of  3 * 10  cm  , and 
—16  2 

E/N  = 1.8  x 10  V cm  . A y.  of  .09  has  been  chosen  and  the  effect  of 

+ 

photons  has  been  neglected  in  the  calculations  shown. 

Figure  6 shows  that  the  electric  field  is  reduced  by  charge  separation 

in  the  region  of  maximum  initial  ionization  at  the  early  time.  The  density 

difference  between  electrons  and  ions  is  very  small  in  the  plasma  region. 

The  field  near  cathode  increases  as  the  electrons  are  swept  toward  the 

anode.  The  field  increases  near  the  cathode  until  it  becomes  high  enough  to 

cause  an  ionization  wave  to  move  toward  the  cathode.  This  is  illustrated 

in  Fig.  7 for  two  different  times.  Note  that  the  figures  show  only  the 

small  region  within  .03  cm  of  the  cathode.’  The  electric  field  near  cathode 

becomes  high  and  forms  the  cathode  fall.  At  t = 0.13  psec,  as  shown  in 

2 

Fig.  8,  the  current  density  is  10  amperes/cm  and  the  thickness  of  cathode 

-3 

fall  is  about  2 x 10  cm. 

The  sequence  of  electric  field  and  charge  distributions  in  the  cathode 
fall  region  is  essentially  independent  of  gap  length  with  the  time  scale 
decreasing  from  ~200  ns  for  a 3 mm  gap  and  100%  overvoltage  to  ~60  ns  for  a 
7 cm  gap  and  50%  overvoltage. 

Figure  9 shows  the  growth  of  the  electron  and  ion  current  density 
through  the  gap  as  calculated  using  the  relation 
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Fig.  5.  Initial  electron  and  ion  densities  and  electric  field. 
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Fig.  6.  Electron  and  ion  densities  and  electric  field  at  0.02  psec  after 

initiation  of  discharge,  showing  depletion  of  electrons  near  cathode 
by  flow  toward  anode. 


DISTANCE  FROM  CATHODE 
(cm) 

Fig.  8.  Electron  and  positive  ion  densities  and  electric  field  near  cathode 
at  0.13  psec  after  initiation  of  discharge.  These  profiles  appear 
to  have  reached  a quasi-steady-state  form.  Note  that  the  calculations 
in  Figs.  5-8  were  carried  out  without  delayed  onset  of  ionization. 

These  curves  are  very  similar  to  those  obtained  with  the  delay 
except  that  in  the  present  case  the  sheath  is  thinner  and  the  voltage 
drop  is  smaller. 
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It  will  be  noticed  that  Eq.  (7)  relates  the  energy  input  to  the  gap 
to  the  integral  of  the  energy  dissipated  and  stored  per  unit  volume 
and  includes  both  the  charged  particle  and  displacement  currents. 

C.  Discussion 

Although  the  mean  free  path  for  electrons  with  an  average  energy  of 
~5  eV  is  small  compared  to  the  cathode  sheath,  a number  of  authors^? 
have  questioned  the  applicability  of  the  continuum  model  implicit  in 


Eqs.  (1)  and  (2).  As  a partial  answer  to  this  question  we  calculate  the 
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energy  relaxation  distance  for  conditions  corresponding  to  those  of  Fig.  8. 


-5 


The  energy  relaxation  distance  for  the  plasma  region  is  about  4 x 10  cm 


which  is  very  small  compared  to  the  dimension  of  the  plasma  region.  A 

,-5 


relaxation  distance  of  3 x 10  cm  is  found  at  the  midpoint  of  the  cathode 
fall  region.  Since  the  electric  field  varies  by  only  4%  over  this  distance, 
we  believe  that  the  continuity  equations  are  reasonably  accurate.  We  have 


compared  our  results  to  the  experimental  values  for  the  steady  state  for 
,29... 


cathode  fall  in  Fig. 10.  We  see  that  the  calculated  cathode  fall  voltage  drop 
is  close  to  steady-state  values  measured  at  much  lower  gas  densities.  The 
calculated  product  of  the  pressure  and  the  cathode  fall  thickness  is  larger 
than  the  low  pressure  values  by  about  a factor  of  two. 

Obviously,  further  work  is  necessary  to  clarify  the  relationships  among 
the  various  theoretical  approaches  to  the  cathode  fall.  At  present,  both  the 
continuum  model  and  the  nearly  free  fall  approaches  seem  to  be  able  to  use 
their  methods  to  show  their  validity.  In  addition,  other  processes  should  be 


included  in  the  model,  e.g.,  the  destruction  and  ionization  of  metastable  atoms 

30 


and  molecules  and  the  heating  and  subsequent  motion  of  the  neutral  gas.' 


4 A 


Fig. 


10.  Exper j mental  and  theoretical  values  of  steady-stale  voltage  drop 

across  the  cathode,  fail.  The  solid  pojnt  Is  the  result  of  our 


calculation!,  using  the  delayed  onset  of  Ionization.  The  curve 
serves  to  connect  i he  experimental  and  theoretical  data. 


-27- 


IV.  ELECTRON  EXCITATION  OF  MOLECULAR  METASTABLES 


Drs.  S.  A.  Lawton  (to  9/77),  D.  Levron  (from  9/77),  and  A.  V.  Phelps. 


The  objective  of  this  project  is  the  measurement  of  excitation  and 


destruction  rate  coefficients  for  metastable  states  of  molecules  which 


serve  as  upper  states  of  lasers  or  which  make  possible  efficient  excitation 


of  such  laser  states.  The  analyses  of  measurements  of  rate  coefficients  for 


the  electron  excitation  of  Oj  molecules  to  the  b £ state  have  been  completed 


and  are  being  prepared  for  publciation.  As  indicated  in  the  previous 


Semiannual  Report  the  measured  excitation  coefficients  for  inter- 


mediate electron  energies,  i.e.,  average  energies  near  2 eV,  are  a factor  of 


three  larger  than  expected  on  the  basis  of  electron  beam  experiments  and 


molecular  potential  data.  If,  as  appears  to  be  the  case,  this  excitation  of 


the  b molecules  is  the  result  of  cascading  from  the  "6  eV  levels' 


then  the  energy  efficiency  for  the  excitation  of  these  levels  must  be  near  60% 


“16  2 

for  E/N  of  about  1.5x10  V-cm  . Recommendations  for  resolving  this  discrepancy 


and  providing  a good  energy  balance  for  electrons  in  include  measurement 


of  the  rate  coefficients  for  vibrational  excitation,  for  excitation  of  the 


O2  Ag  state,  and  for  dissociation  of  C^. 


During  this  report  period  we  have  initiated  measurements  of  rate  coefficients 


for  the  excitation  of  the  A(  E ) state  of  Nj.  This  state  is  well  known  for  its 


energy  storage  properties  in  pure  and  in  N2~rare  gas  mixtures  under  electron- 


beam  and  discharge  excitation.  Problems  being  addressed  in  connection  with 


these  measurements  include  improvements  in  the  cleanliness  of  the  vacuum  system. 


a low  noise  monochromator-photon  detection  system  to  allow  a measurement  of 


the  distribution  of  emitted  radiation  in  the  various  ultraviolet  and  visible  bands 


of  N?,  a reduction  in  the  intensity  of  scattered  radiation  from  the  uv  source 


used  to  provide  the  photoelectrons  used  in  the  drift  tube,  and  the  absolute 


calibration  of  the  photon  detection  system  at  the  ultraviolet  wavelengths  emitted 


3 + 

by  the  A(  E^)  metastables,  i.e.,  the  Vegard-Kaplan  bands. 


t A * 


SCATTERING  AND  TRANSPORT  OF  RESONANCE  RADIATION 


Drs.  T.  Fujimoto  (to  9/77),  A.  Zajonc  (from  9/77),  and  A.  V.  Phelps. 

The  objective  of  this  project  is  to  determine  the  rates  for  both  radiative 
and  non-radiative  energy  loss  from  resonance  states  of  a typical  metal  vapor 
atom  (Na).  The  data  resulting  from  such  measurements  are  necessary  for 
accurate  modeling  of  proposed  high  power  metal  vapor-rare  gas  excimer  lasers. 

During  this  report  period  we  have  completed  measurements  and  most  of 
the  analysis  of  the  spectral  intensity  of  scattered  resonance  radiation 
when  Na  vapor  is  excited  with  white  light.  With  this  mode  of  excitation 
the  incident  radiation  excites  atoms  throughout  the  cell  and  the  dominant 
loss  processes  are  radiative  transport  and  quenching  by  Na2  and  by  any 
added  gas,  i.e.,  by  Nj  in  some  of  these  experiments.  Representative  spectral 
intensity  curves  are  shown  for  pure  Na  in  Fig.  11,  where  the  upper  curve 
is  the  forward  scattering  and  the  lower  curve  is  backward  scattering. 

The  points  are  measured  values  and  the  solid  curves  of  Fig.  11  are  the  results 
of  calculations  using  theory  developed  previously  [D.  G.  Hummer  and 
P.  B.  Kunasz,  J.  Quant.  Spectrosc.  and  Radiat.  Transfer  1^,  77  (1976)] 
for  the  interpretation  of  these  experiments  and  include  the  smearing  effects 
of  the  measured  monochromator  profile.  One  notes  that  the  agreement  between 
theory  and  experiment  is  very  good  except  at  line  center  for  the  forward 
scattering.  The  origin  of  this  discrepancy  is  unknown,  although  it  could  be 
an  underestimate  of  the  non-radiative  transport  of  excitation  to  the  entrance 
window. 

In  order  to  obtain  a quantitative  measure  of  the  non-radiative  contribution 
to  the  loss  Of  resonance  atoms,  we  have  initiated  measurements  of  the 
scattering  of  resonance  radiation  under  conditions  designed  to  emphasize  the 
effects  of  non-radiative  transport.  The  new  feature  is  the  use  of  monochromatic 
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radiation  rather  than  the  white  light  source  used  previously.  Since  the 
absorption  coefficient  is  very  large  at  line  center  this  method  of  excitation 
enables  one  to  establish  a large  gradient  in  the  excited  atom  density. 

Diffusive  effects  are  expected  to  be  important  under  such  conditions. 

To  this  end  a single-mode  argon-ion  pumped  tunable  CW  dye  laser  has  been 
made  ready  for  use. 

Previous  measurements  on  the  radiative  and  non-radiative  transport  of 

resonance  radiation  have  been  limited  to  temperatures  below  500  C 
16  ~3 

( [ Na ] £ 5x10  cm  ).  Since  non-radiative  transport  becomes  more  important 

at  high  vapor  densities,  work  is  presently  underway  to  extend  these  measure- 

18  —3 

ments  to  densities  as  high  as  10  cm  . Because  of  the  difficulties  encountered 
previously  with  "brazed"  sapphire-metal  seals  and  a welded  cell  construction, 
we  are  attempting  to  develop  pressure  seals  with  a thin  metal  gasket 
between  the  sapphire  window  and  the  cell  body.  This  type  of  construction, 
of  course,  requires  careful  design  and  construction  of  the  cell  components 
with  regard  to  coefficients  of  expansion  and  surface  finish.  A prototype 
sodium  cell  is  presently  undergoing  testing. 
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